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Abstract 



In this paper, we made a statistical approach on the thermodynamic structure of 
the nucleon and its quark distributions. We assume that the nucleon is a thermo- 
dynamic system of quarks and gluons. After we derived the quark density of states 
from the Dirac Equation, we calculated each quark's chemical potential, associating 
it with the quark's spin distribution. By using the sum rule of the longitudinal mo- 
mentum, we obtain a function showing the relation between the temperature, radius 
and mass of the nucleon. The nucleon's radius is determined by the total energy 
minimal principle, and thus the temperature is also known. Finally, we deduced the 
quark distribution functions near x=l, and the results are in agreement with the 
experimental data. 



1 Introduction 

It has been justified by several authors in the literature that the nucleon, 
which confines many sea quarks and gluons, could be considered as a thermo- 
dynamic system with temperature T in a spherical volume with the radius R 
[1], [2], [3], [4], [5]. The statistical method based on the grand canonical en- 
semble formalism is then introduced to study the thermodynamic properties 
of the nucleon [6]. After R. S. Bhalerao found a relation between the quark 
distibution functions q(x) in the infinite momentum frame and its density of 
states p(E) in the nucleon rest frame [3], 
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where M is the nucleon mass, g is the quark degeneracy factor and /(//) is the 
Fermi distribution function, the statistical methods can be used to calculate 
various distribution and structure functions of the nucleon [7]. 
However, in order to settle down all the parameters in eq (1), we have to 
find out the quark density of states and their chemical potentials at first. In 
this paper, we deduce the density of states by Dirac Equation and obtain 
the chemical potentials by inputting Au, Ad, As. Our method differs from 
Bhalerao's model [7] in that: 

1. We use a new density of states derived from D.E.; 

2. We determine the nucleon's radius by the T-R-M function and the total 
energy minimal principle; 

3. We use analytic methods, and our results are only about the x ~ 1 region. 
Futher results may able to be extented to small x region if we use numerical 
methods. 



2 Dirac Equation and Quark Density of States 



In our recently proposed Quark Gluon Coupling Model [6], we describe the 
nucleon as confining the quarks and gluons in a spherical volume with the 
radius of R. The quarks interact through the exchange of gluons, and the 
gluon couples to the conserved quark-current. The Lagrangian density for this 
field is 

L = i> \r (id, - g c A») ~ m q ] i> - -F, V F^ , (2) 



where A, is the gluon vector field, g c is the coupling Constance, m q is the quark 
mass and = d,A u — d u A,. If we replace the vector field operators by their 
expectation values, that is, A, — > S^Ao, the Lagrange's equations yield the 
the linear Dirac Equation, 

i-fdp - g c j°A - m,U = 0. (3) 



The normalized quark wave function within a spherical volume with the radius 
of R is 



^(f,t) = Ne x 



( ^ E ~ 9c T mq ^(E-g c A,f-mlr) \ 



/ E-g c A ~m,g 
l\l e 



h{yj {E - gJky - mfr)) ] 



Xq , (4) 
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where r = \ f\, Xq is the quark spinor and N is the normalization constant [8] 
and E is the quark energy eigenvalue. We note that ipt/j should vanish at the 



2 



boundary in the relativistic theory [9], ipip = 0. This yields the energy 
states [6]: 



E n » ^ C A + ^m 2 , + (— ) 



(5) 



where n = 1, 2, and therefore the quark density of states is 



If we use the assumption in the parton model that m q — > 0, the density of 
states will be reduced to p(-E') = — . We will use this density of states in all of 
the following sections, and we find the results are satisfactory. 



3 Chemical Potential and Spin Distribution 



In this section we will discuss the relation between quarks' chemical potentials 
and A-u, Ad, As, R. Many experiments show that we should take quarks with 
spin parallel to the nucleon's spin (q j) different than those with spin anti- 
parallel (q I). The main compoents of the nucleon, except for the gluons, are 
u t, u I, d t, d I, s t, s i and their anti-particles u [, u f, d {, d f, s {, s | 
[7]. We assume that heavy quarks c, b, t almost do not appear [3]. 
In the nucleon, there are 12 constraints on the numbers and chemical potentials 
of the above 12 types of quarks. If q stands for u, d, s quarks generally, Aq 
stands for Au, Ad, As, n q stands for valence quark number and \i q stands for 
chemical potential, we have the constraints as follow: 

(wg + Aq) 

n Ql- n ql = g ' ( ' 

K - A <l) /ON 

n q i ~ ngt = ; (8) 

AM + Vvl = °> ( 9 ) 
Ni + = °- ( 10 ) 



By integrating eq (1) over x, we obtain that: 

n q = / 3f(fi q )p(E)dE (11) 

J 
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Using the density of states deduced in section 2 and assuming that the nu- 
cleon's mass is much biger than /i q , we get: 



3RT 



ln[exp(^) + 1] 



(12) 



n q = 



7T 



From eq (7), (8), (9), (10), (12), we get: 



(rig + Aq)n 
6R 



(13) 



Nl = 



(n q - Aq)ix 
QR 



(14) 



The relation between Aq and chemical potentials is 



Aq 



3R(^ - fi ql ) 



(15) 
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From several experimental data [10], [11], [12], we use Au = 0.83, Ad = —0.43, 
As = -0.10 for the proton and Au = -0.40, Ad = 0.86, As = -0.06 for the 
neutron in the following sections. 

Once the nucleon radius R is determined, all the chemical potentials are known 
from eq (13), (14), (9), (10). Even when R is unknown, we still can get three 
conclusions from those equations: 

1. The ratios of each chemical potentials are: fj, u ^ : fi u i : fi^ '■ f^di '■ /^st : / i s| = 
2.83 : 1.17 : 0.57 : 1.43 : -0.10 : 0.10 for the proton and /z^ : /i u [ : jj,^ '■ Hdi '■ 
^ : ^ sl = 0.60 : 1.40 : 2.86 : 1.14 : -0.06 : 0.06 for the neutron. 

2. For both proton and neutron, u|,uj,,d|,dj,,sj,,sj, have positive chem- 
ical potentials while their anti-particles have negative chemical potentials. 

3. The absolute value of the chemical potentials ofst,sJ,,s~t,sj are much 
smaller than others. 



4 T-R-M Function of the Nucleon 

In this section we will discuss the relation between nucleon's temperature, 
radius and mass and finally determine the value of R and T. 
The longitudinal momentum fraction of each types of quarks could be obtained 
by integrating over x the eq (1) mutiplied x: 




(16) 
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Similarly, the longitudinal momentum fraction of the gluons is: 

4 rM/2 



P9= 3mIo l ^) P {E)EdE, (17) 



where /b(0) is the Bose distribution function and the chemical potential of 
the gluons is 0. Our calculation shows that for both proton and neutron, 
Pui, Pdh Pdi and p g are relatively large, while p sh p si , p s i, Psi are relatively 
small and p u ^, p u ±, p^, p^i are relatively tiny. From eq (16), (17), we have: 

n(n u + Auf 2ttT 2 R 

Pu ' ~ isMR + ~^r ] (18) 

_ n(n u - Auf 2nT 2 R_ 
Pul 18MR + 3M ' 1 ' 

_ n(n d + Adf 2nT 2 R 
Pd ' ~ 18MR + (20) 

= rr(n d - Ad) 2 2nT 2 R 
Pdl 18MR 3M ' [ ' 

32rrT 2 R 

p 9 = (22) 

Psi = Psi = Pn = Pn = (23) 

P^ = Pui = P^ = Pdi ~ 0, (24) 



where n u and n d are valence quark numbers. Note that the sum of all longi- 
tudinal momentum fractions is 1, that is, 

Pul + Pul + Pdl + Pdl +Pg+ Psi + Psi + Psi + Psi + Pul + Pul + Pdl + Pdl = l-(25) 



From eq (18), (19), (20), (21), (22), (23), (24), (25) we get the T-R-M function 
as follows: 

w 68ttT / 5 + Aw 2 + Ad 2 \ . . 

M = RT + . 26 



The total energy of the sphere Et includes the zero point Casimir energy ^, 
where C ~ 2: 

E T = M + ^ (27) 
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We determine the nucleoli radius R by the total energy minimal principle: 

(dE T \ 



OR ) T 



0, (28) 



and therefore get: 



2,(5 + A„3 + A^ + j) 

9M v ; 

9M,/(5 + Am 2 + Ad 2 + 2£) 

47T V / T7(5 + A M 2 + Arf 2 + f )' 1 ' 



For the proton, if we use M=938 MeV, Am = 0.83, Ad = -0.43, C = 2, we get 
i? = {IhAMeVy 1 = 1.28 fm, T = 63.5 MeV. And for the neutron, if we use 
M=939 MeV, Am = -0.40, Ad = 0.86, C = 2, we get R = {IhZMeVy 1 = 1.29 
fm, T = 63.4 MeV. All chemical potentials can be obtained from eq (13), (14), 
(9), (10). Both the proton and neutron have similar radius and temperature, 
although the proton has a slightly smaller radius and higher temperature. We 
may take their radius and temperature the same in some approximations. 



5 Quark Distribution Functions Near x=l 



After we get the R, T and chemical potentials of each type of the quarks, from 
eq (1), we can obtain the quark distribution functions near x=l. We assume 
that x is so near 1 that 1 — x « ^ 0.135. From eq (1), we have 

q(x) = i RM { l- x) U,-M^\ (31) 



So that, 



u(x) = u^(x) + u^(x) + u^(x) + u^(x) 

Au(x) = Uf(x) — u^(x) + Uf(x) — u^(x) 

12RM(l-x) f-Mx\ . n u n . Auix 

exp — — — sinh — — — smn (33) 

Tix 1 \ 2T J 6RT 6RT' v ; 

d(x) = df(x) + d^(x) + d^(x) + d^(x) 
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12RM(l-x) f-Mx\ . n d n . Adir 

exp — — — cosh — — cosh — — ; (34) 

Tix 1 \ 2T ) 6RT 6RT' v ; 

Ad(x) = d|.(x) — d^(x) + d^(x) — d^(x) 

12RM(l-x) f-Mx\ , n d ir , Adn . . 

Tx exp \-*r) smh mr smh 6KT> (35) 

s(x) = s^(x) + S|(x) + s^.(x) + S|(x) 

As(x) = s^(x) — s^(x) + sJx) — s^(x) = (37) 



Thus, from eq (32), (33), (34), (35), (36), (37), we can calculate 

As(x) d(xi s(xi w '] ieQ x p or the proton, if we use the R, T in section 

4, we have ^ffi = 0.774, = —0.425 ( quite in agreement with the 

quark-diquark model result [13], [14] ), ^ = 0, g = 0.216, £g = 0.098, 
when x — > 1. And for the Neuton, if we use the i?, T in section 4, we have 
M£) = -0.397, ^2 = 0.786, ^ = 0, ^4 = 4.83, 44 = 0.466, when 
x -> 1. So that, when x -> 1, F 2 n /F| = 2.155w™(x)/w p (:r) = 0.453, which 
is quite in agreement with the recent experimental data [15] and the pQCD 
result [16], [17]. 



6 Conclusion 



The analytic statistical method in this paper is limitted in the large x region. 
However, the advantage of this method is that it shows the relation between 
the quarks' chemical potentials and their spin distibutions, as well as the 
relation between the nucleon 's temperture, radius and mass, the so called T- 
R-M function. There are three input parameters, Am, Ad, As, which cannot be 
determined without any dynamical approach beyond the statistical mechanics. 
Therefore the statistical method in this paper is only a connection between 
the phenomena, but it is helpful for our understanding. 
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